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In this paper we generalise the main result of a recent work by J. L. Cardy and the present 
authors concerning the bi-partite entanglement entropy between a connected region and its 
complement. There the expression of the leading order correction to saturation in the large 
distance regime was obtained for integrable quantum field theories possessing diagonal scattering 
matrices. It was observed to depend only on the mass spectrum of the model and not on the 
specific structure of the diagonal scattering matrix. Here we extend that result to integrable 
models with backscattering (i.e. with non-diagonal scattering matrices). We use again the replica 
method, which connects the entanglement entropy to partition functions on Riemann surfaces 
with two branch points. Our main conclusion is that the mentioned infrared correction takes 
exactly the same form for theories with and without backscattering. In order to give further 
support to this result, we provide a detailed analysis in the sine-Gordon model in the coupling 
regime in which no bound states (breathers) occur. As a consequence, we obtain the leading 
correction to the sine-Gordon partition function on a Riemann surface in the large distance 
regime. Observations are made concerning the limit of large number of sheets. 
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1 Introduction 



A quantity of current interest in quantum models with many local degrees of freedom is the bi- 
partite entanglement entropy pQ. It is a measure of quantum entanglement between the degrees 
of freedom of two regions, A and its complement, in the ground state |gs) of the model. Other 
measures of entanglement also exist, see e.g. PP-[S]- Such measures occur in the context of 
quantum computing, for instance. Since entanglement is a fundamental property of quantum 
systems, a measure of entanglement gives a good description of the quantum nature of a ground 
state, perhaps more so than correlation functions. For the formal definition of the entanglement 
entropy, consider the Hilbert space as a tensor product of local spaces associated to the sites of a 
quantum system. This can be written as a tensor product of the two quantum spaces associated 
to the regions A and its complement: TC = A A. Then the entanglement entropy is the von 
Neumann entropy of the reduced density matrix associated to A: 

Sa = Tr^/j^log/9^ , = Tr_4|gs)(gs| . (1.1) 

In this work we will be interested in 1-dimensional quantum systems. The entanglement of 
quantum spin chains has been extensively studied in the literature [6]-[Il]. The scaling limit, 
describing the universal part of the quantum chain behaviour near a quantum critical point, 
is a quantum field theory (QFT) model (which we will assume throughout to possess Poincare 
invariance). The scaling limit is obtained by approaching the critical point while letting the 
length of the region A go to infinity in a fixed proportion with the correlation length. It is 
known since some time |15]-|18] that the bi-partite entanglement entropy can be re- written in 
terms of more geometric quantities in this limit, using a "replica trick". It is related to the 
partition function Z„(xi, X2) of the (euclidean) QFT model on a Riemann surface Ain,xi,x2 with 
two branch points, at the points xi and X2, and n sheets cyclicly connected. The position of the 
branch points correspond to the end-points of the region A in the scaling limit. The relation is 
based on the simple re- writing Sa = — hnin^i ^Tr_4/?^, which gives: 

d Zn{xi,X2) , . 

Sa = - hni — . 1.2) 

n^i an z" 

Naturally, this expression implies that we must analytically continue the quantity Zn{xi,X2) 
from n E N, where it is naturally associated to Riemann surfaces, to n G [l,oo). The object 
Tr_4/?^ certainly has a well-defined meaning for any n such that Re(n) > 0. Indeed, Pa is 
hermitian (and has non-negative eigenvalues summing to 1), so that Tr_4p^ is the sum of the n^^ 
powers of its eigenvalues (with multiplicities). Note that this is an analytic continuation from 
positive integers n to complex n that satisfies the requirements of Carlson's theorem p[9j, hence 
the unique one that does. The scaling limit of this object is what defines the proper analytic 
continuation of Zn{xi, X2)- It is natural to assume, as it has been done before [15], that the two 
branch points just become conical singularities with angle 27rn, the rest of the space being flat. 
This is the point of view that we will take. 

In [20], the ratio of partition functions Zn{xi, X2)/Zl'' was studied at large distances X2I = 
r for certain 1-|- 1-dimensional integrable QFTs. The main feature of these models is that there 
is no particle production in any scattering process and the scattering (S) matrix factorizes into 
products of 2-particle S-matrices which can be calculated exactly (for reviews see e.g. [2T]-[25j). 
Taking the S-matrix as input it is then possible to compute the matrix elements of local operators 
(also called form factors). This is done by solving a set of consistency equations \26\ I27j. This 
is known as the form factor bootstrap program for integrable QFTs. In |20j . this program was 
used and generalised in order to compute Zn{xi, X2)/Z^ in the case of integrable models with 
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diagonal scattering matrix (that is, without backscattering) . It was deduced that for this class 
of models, the entanglement entropy behaves at large length r of the region A as 



Here, are the masses of the i particles in the QFT model, with mi < nia- The first term 
is the expected saturation of the entanglement entropy (with c the ultraviolet central charge), 
the precise value of which depends on the details of the microscopic theory included into the 
non-universal small distance e. The interesting feature is the universal second term, where we 
see that the leading exponential corrections are independent of the scattering matrix, and only 
depend on the particle spectrum of the model. This is quite striking: for instance, a model 
of £ free particles of masses will give the same leading exponential corrections as one with 
diagonally interacting particles with the same masses. 

The purpose of this paper is to point out that this result still holds for non-diagonal integrable 
models, and to analyse the particular case of the sine-Gordon model. We obtain the first 
large-distance correction to Zn{xi,X2)/Z^ for that model and provide a detailed analysis of its 
large-n behaviour. This analysis gives a numerical confirmation of the main result ()1.3p . The 
large-n behaviour also shows interesting features that are in relation with the properties of the 
perturbing field: it is linear oc n in the super-renormalisable case, and has an extra logarithmic 
factor oc n log n in the marginally renormalisable case (where the sine-Gordon model represents 
a sector of the 5'C/(2)-Thirring model). 

The paper is organised as follows: In section 2 we review the relationship, discussed in [20] . 
between the partition function (hence the entanglement entropy) and the two-point function of 
branch-point twist fields in a theory consisting of n copies of a given integrable model. We gen- 
eralise the procedure introduced in [20j for computing the entanglement entropy from analytic 
properties of form factors of branch-point twist fields to integrable models including backscat- 
tering. This is based on expressing the two-point function in terms of a form factor expansion. 
For all such models we conclude that, under certain assumptions, the leading order correction 
to the entropy in the infrared limit is of the same form as that obtained in [20] for theories 
with diagonal scattering. In section 3 we introduce the sine-Gordon model and obtain the two- 
particle form factors of the branch-point twist fields. In section 4 we check these form factors 
for consistency by extracting the underlying conformal dimension of the twist fields. In section 
5 we check our general result from section 2 for the sine-Gordon model, and provide a large-n 
analysis of the two-particle contribution to the two-point function. We find good analytical and 
numerical agreement with the results of section 2. In section 6 we present our conclusions and 
outlook. 

2 Form factor equations and the entanglement entropy 

2.1 Form factor equations 

This section follows closely [201. Let us consider some model of 1+1-dimensional QFT. For 
clarity, we will characterise it by its lagrangian density C{(j)) depending on some "fundamental" 
field (j), although the results hold also when no lagrangian formulation is available. In [20], it 
was shown that the partition function Zn(xi,X2) on Mn,xi,x2 is proportional to a two-point 
correlation function of certain twist fields in an extended model, composed of n independent 
copies of the initial theory. The lagrangian density of this extended model, for instance, is 
C^"'\4'i, . . . ,4>n) = YlJ=i This model possesses a natural Z„ internal symmetry under 
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cyclic exchange of the copies, a : = l)---!'^ — l)i (kn ^ 4>i- By a standard 

procedure (see, for instance, the explanations in [20]), to this symmetry and to that of the 
opposite cyclic exchange (T~^, we can associate, respectively, the twist fields T and T, called 
branch-point twist fields. Their two-point function is simply related to the partition function of 
the original model on the Riemann surface: 



Z„{xi, X2) _ ^ ^2d 



^„e2rf"(T(xi)t(x2)) (2.1) 



"1 



where here an below we use (• • •) for denoting correlation functions in the extended model 
Here e is some short-distance scale, Zn is a non-universal normalisation (with Zi = 1), and d„ 
is the scaling dimension of the twist fields [17^ [20] 

dn = ^(n-l], (2.2) 



12 V 

where c is the central charge of the ultra-violet conformal field theory associated to C (that 
describes the short-distance behaviours of correlation functions in the model C). The funda- 
mental property of the twist fields T and T, as operators on the Hilbert space of is the 
"semi-locality" property with respect to any local field O: 

0{y)T{x) = T{x){aO){y) x^ < 

0{y)T{x) = T{x)0{y) x^ > y' 

0{y)f{x) = f{x){a-^0){y) x^ < y' 

0{y)f{x) = f{x)0{y) x' > yK (2.3) 

This implies that O and T (or T) are mutually locafl only when O is invariant under a. Along 
with the fact that they have the lowest dimension, given by (12. 2p . and that they are invariant 
under all other symmetry transformations of the model vC^"^ that are in agreement with a, 
this fixes T and T uniquely up to normalisation (for definiteness, we will assume the CFT 
normalisation, (T(x)T(O)) ~ \x\~'^^"). 

An important remark is that these twist fields are local fields: they commute with the energy 
density at space-like distance^. This is a consequence of the fact that they are associated to a 
symmetry. In particular, the resulting two-point function, which is proportional to the partition 
function of with a defect line extending from xi to X2 through which the fields are affected 
by the symmetry transformation, is independent of the shape of this defect line (this is sometimes 
called a topological defect). This is simply related to the fact that the partition function of C 
on the Riemann surface Mn,xi,x2 is independent from the shape of the branch connections. 

Let us now consider £ to be a massive integrable QFT model. The Hilbert space of massive 
QFT is described by physical particles, and there are two bases, one corresponding to well- 
defined separated particles coming from far in the past, the other corresponding to those leaving 
far in the future. The overlap between these bases is the scattering matrix, and the main 
characteristic of integrable models is that their scattering preserves the number of particles 
involved and the set of momenta, and that it factorises into two-particle processes, as if particles 
were scattering by pairs at very distant space-time points from each other. Hence, only the two- 
particle to two-particle scattering matrix is relevant, and consistency under all possible pair-wise 



*That is, they commute, hence are quantum mechanically independent, at space-like distances. 
^The set of all local fields that constitute a QFT is the set of those that are mutually local with respect to the 
energy density. 
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separations gives the Yang-Baxter equation. Consider an integrable model with mass spectrum 
{rua, a = !,...,£} (some masses may be equal). Let us denote by \0i, ■ ■ ■ ,Ok)ai...ak with 
9i > ■ ■ ■ > 9k the in-states with k particles of rapidities 9i, . . . ,0^ and, respectively, of particle 
types (quantum numbers) ai,...,ak- For 6i < . . . < 9k, this will represent the out-states. 
Then, the two-particle scattering matrix is defined by 

1^1,^2)0102 = X] Sala2i(^l - ^2)1^2, 6'l)/32/3i 

/3i,/32 

where we used relativistic invariance in order to write the scattering matrix as function of the 
rapidity differences. This equation holds for 9i > ^2 as well as 9i < 92, thanks to unitary. One of 
the main achievements of massive integrable QFT is the exact calculation of scattering matrices 
in many models, from its expected physical properties along with the Yang-Baxter equation (or, 
in some cases, from Bethe ansatz solution of an underlying integrable microscopic model). 

Another important result of massive integrable QFT is the exact calculation of some form 
factors in many models [26\ [27| (see [28] for a recent review) . Form factors are matrix elements 
of local fields between the vacuum and a many-particle state (say an in-state): 

F^^'''-'"'{9i, . . . ,9k) = {vac\Om9i, . . . ,9k)a^...a, • (2.4) 

They are, more precisely, the analytic continuation of these matrix elements in the rapidity 
variables. The exact evaluation of these objects follows from solving a set of expected properties 
that involve the exact two-particle scattering matrix and form a Riemann-Hilbert problem, 
along with certain minimality assumptions on the analytic structure. Besides integrability, 
an important requirement for this Riemann-Hilbert problem to hold is locality of the field O. 
Indeed, a strong indication of its validity is that it is possible to show that two fields whose form 
factors solve it commute at space-like distances. For the simplest local fields, this Riemann- 
Hilbert problem is well known and can be solved in many cases. 

If C is integrable, then certainly also is, with n times as many particles, which we will 
denote by the double index (a, i) for a = 1, . . . , £ and i = 1, . . . ,n. Its scattering matrix is 
simply given by 




h / «2 
k = i2 



Locality of the twist fields T and T along with the exchange relations (12. 3p were used in |20j in 
order to justify a Riemann-Hilbert problem for their form factors in integrable models with diago- 
nal scattering (that is, with Salal{9) = 80102(0)^01,13-1^^02,132) ■ The generalisation to non-diagonal 
scattering is straightforward. We will consider only T, since T = T'^ on the Hilbert space. Let 
us denote by /x, v and uj double indices of the form (a,j), and by f'^^^^'"^'' {9i, . . . ,9k, n) the 
form factors of T in the n-copy model . They are analytic in the rapidity differences except 
for poles (that may accumulate at infinity), and satisfy the relations 

i^k i- • • ,Pi,t/j+i, . . . ,nj — b^.fj,^^^{Uii+i)l<,^ {. . . . . . ,n), 

1^1,1^1+1 

^rK;.2.../^..(^^^27ri,...,0fc,n) = F^''^-''^''' {92, . . . ,9k,9i,n), (2.6) 
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and 



-,T|i/i...i/fc 

with 



j/i,...,j/fe 

xFP'''-'"'{9,,...,ek,n) (2.7) 



Here = 9i—9j. As function of 9i for real 92, ■ ■ ■ , 9k, there are no poles of F^''^^ "'^* (^i, . . . , 9k, n) 
in the strip lm{9i) G [0,7r], except for those given by the last two equations, and for poles at 
purely imaginary values (and with purely imaginary residues) corresponding to bound states 
between the associated particles (note that we will not consider any bound states in the sine- 
Gordon example studied below). In the second equation, the crossing or locality relation, we 
introduced the symbol fi = {a,j + 1). As compared to the usual form factor equations, it is 
altered by the nature of the exchange relation and it now relates form factors associated to 
different particle sets (belonging to different copies). Finally, the last two equations generalise 
the standard kinematic residue equation to branch-point twist fields, where we introduced the 
symbol p, = (a, j) with a denoting the quantum number of the anti-particle of a in the theory 
C. Once more, the exchange relations (j2.3p are responsible for the splitting into two equations, 
and the shift in fi. 

It is instructive to specialise to two particles; this is what gives the main result (jl.3p for the 
entropy (one-particle form factors of spinless fields are ^-independent). In this case, the first 
two form factor equations specialise to 

pT\i»,j)im^O^^^ = J]F2^l(^''=)(^'^-)(-0,n)5;^;f);j;t)(^) = FP^^''^^^''+'\27r^ - 9,n), (2.9) 

(5,7 

for all values of j, k, a and (3, where 9 is now the rapidity difference. From the equations above, 
from (j2.5p and from application of the Z„ symmetry, it follows: 

pT\{a,i)((3,i+k)^^^^^ = F2^l("'^')(^'^'+'^)(0,n) V i,j,k,a,f3 (2.10) 

^T|Ki)(/3j)^^^^^ = F2^l('''')("'')(27r(j-l)i-0,n) V a, j e {2, . . . , n + 1}.(2.11) 

The last equation at j = n + 1, and the first equation of ()2.6p . give 

^r|(a,i)(Ai)^^^^^ = ^F2'l(^'i)(^'i)(-0,n)55(0) = F2'l(^'^)("'^)(27rm - 9,n), (2.12) 

(5,7 

so, as in [20], we can solve for the form factors associated to particles in the first copy and then 
use (12. lip to obtain all other solutions. From now on we will abbreviate 

F^\^^^'^''^^'\9,n):=Fj\'^\9,n). (2.13) 

The simple form of the scattering matrix for particles living in different sheets, and the fact 
that no bound state can occur between particles in different copies, simplifies drastically the pole 
structure of the form factors. A combination of the equations above along with the two-particle 
case of the last two equations of ()2.7p gives that f'^^"'^ {9,n) is analytic for lm(0) G [0,27rn] 
except if a = /3 for poles at = ivr, with residue i{T), and at 9 = 2i7rn — iTT, with residue —i{T), 
and for possible poles in lm(9) £ (0, vr) and in lm{9) £ (27to — vr, 27rn) corresponding to bound 
states. 
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2.2 Conical singularities and the entanglement entropy 

The analytic structure described above is sufficient to establish the result ()1.3p . following 
arguments of [20j. The most delicate point of these arguments is the analytic continuation in 
n, which still needs further justification. Below we attempt to support the arguments from the 
geometrical picture of conical singularities. Another point is that we must assume that form 
factors at real rapidities vanish faster than (n — 1)2 as n — > 1. They certainly do vanish as 
n — >■ 1, since then the branch-point twist field becomes the identity field. They were observed to 
vanish proportionally to n — 1 in ^D] . We first write the two-point function in the two-particle 
approximation : 

2 [ 1-part. ^r- /■ f ^n^^nf fn^ „N -Km^coshei+m^coshe^) 



(2.14) 



a,l3=l_ 



-00 —00 



where 



n 



{T?LAO,n) = (2.15) 



2 



Fl\-\e,n) +j;|F2^l°/^(-0 + 2vrzj,n) 

There is no contribution from the one-particle terms when we take the derivative with respect 
to n and evaluate it at n = 1, since they are squares of (analytically continued) one-particle 
form factors, which are ^-independent and vanish faster than n — 1 as n — > 1 by assumption. 

As for the function fafii^^n), coming from two-particle form factors, let us denote by 
fa,i3{0T''T') its analytic continuation from n = 1,2,... to n G [l,oo). The analytic continua- 
tion in n of the form factors themselves, for fixed rapidity, is natural from the geometrical 
picture of conical singularities with angle 27rn. For any n real and positive, form factors should 
have kinematic poles at 9 = in and 9 = 2iTTn — in representing particles going in a straight line 
on either sides of the conical singularity. They also possibly have bound state poles as described 
above corresponding to bound states forming on either side, and they satisfy all other properties 
stated above, now with n real and positive. For technical reasons, we must assume that the 
residues at bound state poles are not diverging as n ^ 1. It is natural to expect that for any 
fixed rapidity, the resulting form factors do not "oscillate" as functions of n. More precisely, 
that they have definite convexity, at least for large enough n. As for the sum over j in (j2.15p . it 
should be understood, when multiplied by n like in (|2.14|) . as the total two-particle contribution 
with particles allowed to cover the region all around the conical singularities. Intuitively, it 
should be smooth and should not present any oscillations either in n, since the space around 
the conical singularities just increases linearly with n. Hence it should have definite convexity, 
at least for n large enough (where parallel incoming straight trajectories just to the right and to 
the left of the conical singularity are far apart past it). In fact, a natural and more constraining 
requirement is that no oscillatory terms are present in the large-n expansion. These properties 
may fix uniquely the analytic continuation in n. 

A smooth function of n > 1 can be obtained easily from the sum over j. Methods used in 
[20j lead to the fact that in the limit n ^ 1, the only possible contribution to the derivative 
w.r.t. n is from the collision of the kinematic poles in f'^^"'^ {9 ,n) when a = ]3. This uses the 
assumption that form factors vanish faster than (n — 1)2 , and the fact that these are the only 
singularities that can collide as n — > 1. Hence the contribution can be evaluated exactly by 
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I I 2 1 /" 

^\pr\aP(^_0 ^2TTij,n) = — (b dsg{s) cot ITS- Q (2.16) 



extracting these kinematic poles. For completeness and clarity, let us provide here the explicit 
analytic continuation of the sum over j, following [20] and emphasising the general principles 
used. We write 

2 

2i 

where 

g{s) = Fp'^'^i-e + 2ms,n){F^^'"^)*{-e -27ris,n) . (2.17) 

The contour ranges from Re(s) = e to Re(s) = n — e encircling counter-clockwise the segment 
s £ [e,n — e] for some < e < 1/2 near enough to 1/2 so that no bound state singularity can 
be encircled. If the form factors at infinite rapidities vanish sufficiently fast (in all examples 
calculated, they vanish exponentially, which is sufficient), then the contour can be widened up 
to Im(s) = iboo, since the contribution at Im(s) = ±00 vanishes. The quantity Q is the sum of 
the residues from the kinematic poles in the product of form factors themselves: 

Q = ReSs=s(^Fp'''^{-9 + 27ris,n){F2^'"^)*{-9 -2TTis,n)^7Tcot7rs 

kinematic poles S 



6a,ptan Im (^Fp''"{-2e + in) - Fp''''{-29 + 2mn - m)) 



(2.18) 



Using the n-periodicity of cot vrs, the part of the integration at Re(s) = n — e can be re- written, 
and the contour integral in ()2.16p becomes 

-/ dy {g{n + iy - e) cot TT{iy - e) - g{iy + e) cot TT{iy + e)) (2.19) 
^ J-00 

This expression is analytic in n; as n is varied continuously down to n = 1, no pole crosses the 
y integration line. The expression for Q is analytic in n for n > 1 as well. 

Note that we used the function Trcotvrs, which is a function that remains unchanged under 
s 1-^ s + n for any integer n and that has simple poles with equal residues, in the band Re(s) S 
(0, n), only on the integers 1, 2, . . . , n — 1. The poles were used to reproduce the correct sum, and 
periodicity was used to cure the part of the integration contour at Re(s) = n — e, which otherwise 
would have crossed poles as n is varied continuously. However, we wish to emphasise that no 
periodicity property of the function g{s) (related to the form factors) itself is used, contrary to 
the derivation presented in [20]. The function g{s) does not cause problems since as n is varied, 
its poles move accordingly. Could we replace vrcotvrs by another function u{s)7 If we require 
that u{s) do not depend explicitly on n in addition to the properties above, then it is uniquely 
fixed to TT cot vrs. Relaxing periodicity to, for instance, n(s+l) = v{s)u{s) would impose v{s) = 1 
for s integer (in order to have the correct residues), and conditions on convergence at infinity 
would require v{s) = 1. Note also that the expression (j2.19p does not depend on e even for n 
non- integer. 

The integral ()2.19p have definite convexity as function of n, at least for large enough n. This 
follows from expected convexity properties of the form factors themselves. The same holds for 
the residues Q. Integration over 6 should preserve these properties for large enough n, and it 
seems natural to expect that the full large-n expansion does not contain oscillatory terms. Of 
course, a more complete analysis of the n-dependence of Q and of the integral ()2.19p would be 
useful. 

It is clear that the integral in ()2.19p vanishes when n — > 1 faster then n — 1 by our assumption, 
and that Q vanishes at least like (n — 1)^/^ for any 6^0. The only possible obstruction in Q is 
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when — > 0, due to the kinematic pole of the first term in (j2.18p at argument (2n — l)i7r, and 
the kinematic pole of the second term at argument in. Extracting these poles gives, for n ^ 1 
and 9^0, 

n) = fc,./„,(o. 1) ( .(.vt,;-!)) - 2(;-t(; - 1)) ) 

with 

/„,^(0,1) = ^. (2.21) 

Here and below we adopt the convention that fa,aiO, 1) = lim„^i fa,aiO,n). This immediately 
gives 

n)^ _^ = ^S{d)d^,f, . (2.22) 

This has a non-zero measure at 9 = 0, and proves (II. 3p . 

As we mentioned above, there are two main requirements for the argument to be valid: 
first, an appropriate analytic continuation has to be taken, in agreement with the intuition from 
conical singularities, and second, the analytically continued one- and two-particle form factors 
have to vanish, for any real rapidity, faster than (n — 1)2 as n ^ 1. As we said above, for the 
latter requirement, in the cases constructed in [20j, it was observed that form factors vanish 
like n — 1, and this is what is expected in general. Below we provide further evidence for both 
requirements in the case of the two-particle form factors by constructing them explicitly in the 
sine-Gordon model. In particular, we verify numerically the striking fact that fa,a{0, 1) = ^ 
is in agreement with nfa^a{0,n) having no oscillatory terms in its large-n expansion. We also 
find that it diverges like n as n — > 00 in the super-renormalisable case, and like nlogn in the 
marginally renormalisable case. 



3 The sine-Gordon case 

Let us now consider the specific sine-Gordon case. This model can be defined by the lagrangian 

^ = ^ido^? - lidiipf + ^cosiP^) , (3.1) 
where /U has scaling dimension 2/3^ — 2 and /? does not renormalise. We will use the variable 



(3.2) 



For u > 1, the spectrum of the model is known to be composed of two particles with equal 
masses, which we will label by + and — , representing the quantised version of soliton and anti- 
soliton in the classical theory. For simplicity we will consider this region only. For v < 00, the 
model is super-renormalisable. At = 1, the particles are free with fermionic statistics, and the 
model is equivalent to a massive Dirac theory. The model with v > 1 can be seen as the massive 
Thirring model, a perturbation of the massive Dirac theory that preserves the U{\) symmetry. 
From this viewpoint, the asymptotic particles are the positively and negatively charged version 
of the same particle. Ks v ^ 00, the model specialises to (a sector of) the S'?7(2)-Thirring 
model, which is marginally renormalisable. 
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The scattering of soliton with anti-sohton is non-diagonal. The scattering matrix is given by 



smh(^) ^ sinhfHE) 



sinh 



sinh (^) 



with 



5+1(0) = SZZ{e) = -exp 







t sinh(tz^/2)cosh(t/2) 



(3.3) 



(3.4) 



Since branch-point twist fields have zero U{1) charge, their one-particle form factors vanish. 
On the other hand, the relations (j2.12p give the following system of coupled equations for their 
two-particle form factors: 



,n 



,n 



Fp^\-9,n)S%{9) = Fp-^{27rin - 9, 



7,5=± 



Fp^\-9,n)Sl^_{9) = Fp+-{27rtn - 9, 

7,5=± 



n 



n 



(3.5) 
(3.6) 



^1 ^ 

These equations can be diagonalised by introducing the linear combinations F2 {9,n) =b 

F2 {9,n). However, the branch-point twist fields are also invariant under charge conjugation 
+ ^ — . Indeed, this is a symmetry of the action, which keeps unchanged partition functions 
on Riemann surfaces, and branch-point twist fields are associated to such partition functions. 



Hence Fr 



F. 



'{e,n) 



F. 



9, n) and we are left with 



,n 



■.TV 



{2'Kin 



,n 



(3.7) 



It is possible to find integral representations for the combination of S'-matrix elements above, 



SXz{9) + S^+{9) 



sm 



- exp 



dt 



2 sinh 



cosh 



sinh(f) 



where we used 

sin ^{a + x) 
sin |(a — x) 



exp 



t sinh(i/t) cosh(t/2) 

dt 2 sinh t{\ — a) sinh(ta;) 



sinht 



for < a < 1. 



(3.8) 
(3.9) 

(3.10) 



Following [20], it is now straightforward to find a minimal solution to ()3.7|) (that is analytic for 
Im(^) G [0, 27rn]) up to normalisation: 



,n 



-i sinh 



-) 

2n J 



exp 



00 2 sinh 



2 



cosh 



t(i/-2) 



sm 



2 it 



\ ITT J 



sinh(nt) sinh(ti/) cosh(t/2) 



. (3.11) 



The two particle form factor F2 {d,n) can be fixed by including the right pole structure as 
in [20]: 

(T)sin(^) Ff^{6,n) 



'{9,n) 



^ 2nsinh(^)sinh(^) F2™'^(i7r,n)^ 

where the normalization has been chosen so that the kinematic residue equation gives 



iT 




(3.12) 



(3.13) 
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4 Identifying the ultraviolet conformal dimension 

It is interesting to check the solutions above for consistency. A possible consistency check 
consists of analyzing the short-distance behaviour of two-point functions involving the twist 
field and compare that behaviour to the one expected from conformal field theory predictions. 
In particular, one can look at the two-point function of the twist field with the trace of the 
stress-energy tensor, Q. It was found in [30j that 



A 



r 



A 



r 



1 



2{T) 



r (e(r)T(O)) dr 



(4.1) 



(where the integration is on a space-like ray). This formula is known as the A-sum rule. The first 
equality, expected from CFT, holds from the A-sum rule thanks to the fact that @ commute 
with T and that ©1^ = 0. Following the derivations in |17[ [20] the conformal dimension is 
expected to be A'^ = c?„/2, where d„ was given in (12. 2p . 

As explained in detail in [20] we can now employ the expansion of the two-point function in 
terms of two particle form factors and obtain 



A 



r 



n 



e|+- 



n2 



h2,n) 



— oo — oo 



n 



2(27r)2m2 (cosh 6*1 + cosh 62)^ 

' deF^^^-{e)Fp^-{e,n)* 



cosh2(0/2) 



(4.2) 



where we have used the fact that Fr 



01 



n2 



f: 



a I 



'12 J 



and F^^^^{9i2) = 0. 



The two particle form factors of G can be obtained from results already known in the 
literature. In particular, by employing the results obtained in [311 [32j for the form factors of 
exponential fields and the relation between the sine-Gordon coupling constant and the soliton 
mass we find: 

e|+_^^^ 27r im? G (6) cosh ^ 



'(0) 



where 



G{0) 



-i sinh — exp 



V sinh(HL^) 

°° dtsmht{l - v)sm^ it (l - ^) 
t sinh(2t) cosh t sinh(tz^) 



(4.3) 



(4.4) 



A table of values of the integral (|4.2p for different values of n and v is presented below. The 
values given in the top row in brackets are the exact values of the dimension, as predicted by 
CFT. Recall that we assumed the non-existence of bound states, which is why we only consider 
the case v > \. 





n = 2 (0.0625) 


n = 3 (0.1111) 


n = 4 (0.1563) 


n = 5 (0.2) 


V = 1.1 


0.0627 


0.1115 


0.1568 


0.2008 


u = 1.6 


0.0628 


0.1118 


0.1573 


0.2013 


u = 2.1 


0.0627 


0.1115 


0.1570 


0.2010 


u = 2.6 


0.0626 


0.1113 


0.1566 


0.2003 


u = 3.1 


0.0625 


0.1111 


0.1562 


0.1999 


= 3.6 


0.0624 


0.1111 


0.1560 


0.1995 
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n = 6 (0.2431) 


n = 7 (0.2857) 


n = 8 (0.3281) 


n = 9 (0.3704) 


n = 10 (0.4125) 


V = 1.1 


0.2440 


0.2868 


0.3293 


0.3717 


0.4139 


V = 1.6 


0.2447 


0.2876 


0.3303 


0.3728 


0.4152 


V = 2.1 


0.2443 


0.2871 


0.3297 


0.3722 


0.4145 


v = 2.6 


0.2434 


0.2861 


0.3287 


0.3707 


0.4129 


u = 3.1 


0.2428 


0.2854 


0.3277 


0.3698 


0.4119 


V = 3.6 


0.2424 


0.2848 


0.3270 


0.3690 


0.4109 



As we can see, all values are extremely close to the exact ones, even in the two particle ap- 
proximation. This provides a consistency check for our form factor (|3.12|) . A further check 
concerns the correct normalization of the form factor ()4.3p . Indeed the normalization employed 
here corresponds to -F^'"*" {iir) = 27rm^ and differs from that in other places in the literature 
|31j . A consistency check for this normalization consists of extracting the central charge c = 1 
of the underlying CFT by employing Zamolodchikov's c-theorem [33]: 

3 

'=2 



drr^ {@{r)Q{{))) . 



(4.5) 



Employing again the expansion of the two point function in terms of form factors, changing 
variables and performing one integral as in (|4.2p we obtain 

2 



3 
8^ 



e|+- 



{0) 



-de. 



10 cosh^(0/2) 

This integral can be easily evaluated numerically for different values of u and gives 



(4.6) 





u = 1.1 


u = 2.1 


u = 3.1 


u = 4.1 


u = 5.1 


c 


0.9999 


0.9918 


0.9888 


0.9877 


0.9874 



5 Computation of /±ip(0, n): Numeric and analytic results 



In this section we wish to study the behaviour of the function /±=p(0,n) (see (|2.15p ) in detail 
and, most importantly, show that /±::f:(0, 1) = 1/2 as expected. Let us introduce the notation 

/(n) = /-I (0,n) = / ^(0,n). This function, by definition, is the natural analytic continuation 

of /_| (0,n) from n = 2, 3, . . . to n G [l,oo), that is with the prescription of smoothness and 

monotony as described in paragraph 2.2. We will obtain the form of its full large-n expansion, 
with explicit first few coefficients, showing that it grows proportionally to n for any finite z/, and 
to nlogn for 1^ = 00. This large-n expansion is interesting in its own right, and turns out to 
differ quite dramatically from the 1/n expansion found for the sinh-Gordon model |20| . We will 
then perform a numerical study to confirm it, and most importantly, verify the value /(I) = 1/2. 

Let us commence by our analytical considerations. Note that the first term of (j2.15p is 
zero at 9 = 0. For the summand, the second identity of ()3.7p implies that F2 (2'/rij, n) = 
F2 {2iTi{n — j),n) for j = 1, . . . ,n — 1. In addition, for any fixed value of n the value of 
the function is largest for smaller values of j and decreases quickly as j approaches [n/2]. This 
behaviour is particularly extreme for n very large. In that case, only for j <^ n does the function 
above have non-negligible values. Therefore, the value of f{n) for n — > 00 can be obtained by 



replacing the sum Yl^=i 



'(27rij»p by 



2 V lim 

^ ^ 71. — ¥C\r\ 



\Fr 



(2^ij,n)|2 =/(oo). 



(5.1) 
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This can be written as 



where 



/(oo) 



6(1/ 



a{j) = Hm a{j, n) 



2(2j) 



8 1_ 

2 2u 



with 



and 



a{j, n) 



sin(^)sin(^)^ 



3_ J_ 
1 2 2i/ 



n ■ /(l-2i)7r\ • /(l+2j)7rN - / tt \ ^ ^ 

2nsin( ^ )sm{ ^^2^' )sin(^)2 2. 



(5.2) 



(5.3) 



(5.4) 



b{j) = Hm b{j, n) = hm Fp''{2TTij, n) 



sm 



JTT 

n 



21/ 2 



exp 



00 (2 sinh ( ) cosh ( 



2t 



smh.{tu) cosh(t/2) 



-2tj 



(5.5) 



with 



b{j, n) = exp 



/2sinh(ii^ 



cosh 



t \ sinh(tz^) cosh(t/2) 



J__ 3 

2i/ 2 



u-l] sin^ f (n - 2j) 
sinh(nt) 



(5.6) 



Notice that the factor sin (j^^ ^" ^ in ciU) is cancehed by a similar factor in 6(j). They have 

been introduced in order to guarantee that the integral in b{j) is convergent as t — > 0. Formula 
(|5.2|) holds only for v finite; we will come back to the case = 00 below. Evaluating (|5.2|) 
numerically we obtain the values listed in the table at the end of this section. 

However, contrarily to the sinh-Gordon example studied in [20], the large-n corrections to 
this linear behaviour are not just of 1/n type but depend on the particular value of v. The form 
of these corrections can be determined by some analysis. Notice that 



fin) 



n/2 



6(1/2, n)2^^ 

(«-l)/2 



^ a{j, nfb{j, 



1 



n 



6(l/2,n)2 



6(1/2, n)2 



^ a{j,nfb{j, 



n] 



a{n/2,n)^b{n/2,nf' n even 



n odd 



(5.7) 



Let us first analyse b{j,n) and a{j,n) independently. It is trivial to see that limj^oo6(j)^ = 1. 
More precisely, it has a large-j asymptotic expansion of the form 



buf = Y.^k + o{rn, with 



1. 



(5.8) 



fc=0 



where we use powers of 2j instead of j for later convenience. This formula is obtained from 
a small-i expansion of the coefficient of e"^*-' in the integrand in the last line of ()5.5p . The 
corrections for finite but large n can be obtained by looking at b{j,n)/b{j). Expanding the 
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resulting integrand in t while keeping intact the factor sin^ ^{n — 2j)/sinh(nt) and the factor 
e~^*-^ , then by integrating this expansion term by term, we find 

where (^{z, a) = X]fe^o(^ ~^ a)~^ is the generalized Riemann zeta function. The full series is just 
the exponential of a series that has terms in n~^*^ for k = 1,2,..., with coefficients that are 
functions of j/n, linear in the functions C ^2A;, 1 — + C (^2k, 1 + ^^d- This series expansion is 
convergent for all j = 1, . . . , n/2. 

As for the function a{j), its square has a large-j series expansion 

4 

a{j? = T-T^cikr^'', with ao = l. (5.10) 

7r2(2j) +^ 

The square of its corrected form a{j, n) can be expanded at large n as follows: 

where Pfc(j^) is a polynomial of order k in j^. Again, this series expansion is convergent for all 
j = 1, . . . ,n/2. In fact, it converges for all values of n such that \n\ > \j it 1/2| and \n\ > 1/2. 

We now consider the large-n expansion of ()5.7p . We will do explicitly the case where n is 
even, and comment about the agreement with the case where n is odd afterwards. First, the 
subtraction 0,(^/2, n)^6(n/2, n)^ for n even can be seen, from the analysis above, to be 



-1 



n-'--{l + n-'[[n~% (5.12) 



6(1/2)2 

where here and below we use the notation 

[[x]] = some series in non-negative integer powers of x. (5.13) 

Let us then consider only the term which is a sum, and subtract for convenience the value 
6(1/2)2/6(1/2, n)2/(oo), which is just [[n-2]]. We are left with 

^^^^ a{j, nfbij, nf - ^ a{jfh{j)^ . (5.14) 

The pre-factor 2/6(1/2, n)2 is also so we will consider the sum without pre-factor. We 

will consider three parts: 

n/2 oD oc 

I. Y.a{j)H{jf-Y,a{jfb{jf = - <jfb{j)\ (5.15) 

j=l j=l j=n/2+l 

n/2 

II. ^^{corrections to a(j)2} • 6(j)^, (5.16) 
i=i 

n/2 

III. y^a(j, re)^ • {corrections to 6(j)2}, (5-17) 



*Note that all such functions can be written in terms of trigonometric functions, for instance <^(2, 1 + 2) + 
(^(2, 1 — z) = Tx^ / sin^ TTZ — 1/z^, but the analytic structure is clearer using the generalised zeta functions. 
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where the correction terms are those from (jS.lip and (j5.9p . In the part I, we can use the large-j 
expansions of a(j)^ and b{j)'^. Consider the identity 

oo / 1 1 \ 

E Pi)-- = 1 + f ) = + + n-[|n-ll) (5,18) 

j=n/2+l \ V y / 

(here, the coefficients of the series [[n~^]] are polynomials in z), where Ci^) = Ci^^'^) is the 
Riemann zeta function. This implies that 

I:- Yl <jfKjf = ^[-^ + -+n-'[[n-']]]. (5.19) 

j=n/2+l TT \ n J 

The contributions II can be analysed using the formula 

"/2 .w 00 

E I^JJ _ r,z-2-2fe jj{2+2k-z) 

4p -1 ~ rijl 



j=l •' k=0 

00 



fe=0 

Here, Hj^^ is the Harmonic number of order b 



E ^-'-''({2 + 2k-z)+ n^-i (^(^ + ^ + n-\n-']]^ . (5.20) 



^^''^ = E4 = CW-C(^« + 1) ■ (5-21) 



3=1' 

Note that the infinite sum on the right-hand side of (I5.20p is convergent for any non-integer z. 
Also, the coefficients in [[n~^]] are rational functions of z. Let us write 6(j)^ = hk'{2j)~'^^' + 

r{j) where r(j) is the rest. The function r{j) behaves proportionally to as j — > 00, so 

that the sum from 1 to n/2 of the product of r(j) times the k^^ correction to a(j)^ in (15. lip can 
be extended to a sum from 1 to 00 without problems. Hence for the k^^ contribution to II we 
can write: 

^ Z^Z^bk' Pk{j )+n 2^^0) (^.2_i) PkO ) 

fc'=Oj=l j=i ^ •' ' 

00 CO ,„ .s3_i_2/t' 

-n--- E ^^'%r^Vk{f) . (5.22) 

fc'=fe+l j=n/2+l ^ ' 

The first term can be evaluated using ()5.20p . The second term is a contribution of the order 
^-2k-2 only^ and the last term can be evaluated using the negative of (I5.20p without the infinite 
sum of zeta functions and is of order fi~'^k-2~i/u _ Hence, the result is a term proportional to 
^-2k~2 pj^g g^j^ infinite series of the form n"~[[n^^]]. For instance, the leading contributions to 
this series are 

where pk^k is the coefficient of j"^^ in PkW^)- When we consider all corrections, for all k = 
0, 1, 2, . . ., we obtain, for the coefficient of any given order 7^-1/^+'^ with fixed k' = 0, 1, 2, . . ., 
an infinite sum over k. for instance, these are the sums over all k of (j5.23p in the cases k' = 
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and k' = 1. We must make sure that all these infinite sums give finite results. By putting 
n = aj for a > 1 in (jS.lip . a value where the series in (jS.lip converges for all j > 1, we 
see that YlT=oP>'(j'^)(^j)"'^^ ^ This implies that Y^^=QPk,k-iO'~'^^ < oo, where Pk,k-i is the 
coefficient of j^^-^i -^^ pfe(j^). Since the coefficients in [[n~^]] in (I5.20p are rational functions of z, 
a little analysis of (|5.22p shows that this is sufficient to prove that the infinite sums that are the 
coefficients of 71"^/^+'^ for any k' = 0,1,2, .. . give finite results. For instance, the contributions 
to n~^^'^ and n~^^'^~^ can be found from 



^ n^oo y (2j) ~ 



k=0 



j=ns 



TTS ^ 

sinvrs/ / vr^s^ 



(5.24) 



by 



00 00 r) — 2fc /•- 

y: = Ml/2) , = 2'^' / ' • ^'-'^^ 

fc=0 fc=0 ' -^^ 

The result is that the leading contributions for the case II are 

Finally, a similar analysis can be done with the part III, involving the correction terms of b{j)'^ 
in (j5.9p . by expanding the generalised zeta functions in powers of j/n. Considerations similar 
to those above, with the fact that the expansion of the generalised zeta functions is valid for 
\n\ > lead to the same structure for the large-n expansion as that of parts I and II, with an 
extra factor n^"^^ for each correction term with /c = 1, 2, 3, . . .. This gives: 

III : n-^[[n'^]] + n-^-^[[n-^]] . (5.27) 

Putting everything together, along with the subtraction (|5.12p specific to the case n even, 
we find 

n/(n) = Con + — + + . . . + n"- (don + di + — + . . .] (5.28) 

where 



n n-^ \ n 



Co = /(oo) 



di = . 



It is striking that although the intermediate steps of the analysis give terms n''~'^ for k = 
1,0,— 1,-2,..., the term with A; = identically vanishes for all ly. In fact, one can see that the 
same analysis for the case n odd, starting from the second form of ()5.7p . directly gives vanishing 
coefficients for k = 0, —2, —4, .... For instance, changing n ^ n — 1 into ()5.18p erases the term 
l/(2n) on the right-hand side. The other terms computed above are also unchanged in the odd 
case. Hence we conjecture that the large-n expansion valid from both n even and n odd is 

n/(n) = Con + — + + . . . + n""^ (don + — + ^ + . . . ) (5.30) 

n n'^ \ n n'^ I 
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Indeed a very precise numerical fit of nf^ (0, n) for several values of i^, is given by the 

function above. Let us commence by evaluating nf{n) for integers n > 1 and several values of v 
(see Fig. [1]). Notice that u = 1 corresponds to the free Fermion point and, indeed for this value 



nf .(0,n) 



V = 1 

V = 1.5 

V = 4.2 

V = 9.6 

V = 20.4 

V = infinity 









—t — ' 1 — 

10 20 



~^ — — ^ 

30 n 40 



— I — ' — I — 

50 60 



Figure 1: The function (0, n) computed numerically for several values of v and 2 < n < 70. 

of v we recover the result obtained analytically in [20j for the Ising model, namely /(n) = 1/2 for 
all values of n (in particular, also for n = 1). For other values of i' the large-n linear behaviour 
seems also apparent from Fig. [1] and is clearer the smaller the values of ly. In fact, for = 9.6 
and ly = 20.4, the correct linear behaviour can only be seen at much larger values of n, because 
of the importance of the term in n^~^^'^ . 

Our numerical analysis has revealed that the numerical values of n/-| (0,n) (especially for 

small n) are best fitted by a function of the form: 

F(z/, n) = (1 + n)co + ci + + do{l + n)^-^ + + ^H— ^, (5.31) 

1 + ^ (1 + n)^ (l + n)^+^ 

This function has the same large-n behaviour as ()5.28p but since the expansion ()5.28p is only 
asymptotic its working is not guaranteed for small values of n (in particular at n = 1). Shifting 
n — > 71 + 1 constitutes a re-summation that allows us to consider small values of n. The tables 
below show the values of all constants involved in (j5.3ip as well as the values of cq and do as 
obtained by numerically evaluating the sum ()5.2p and the function (I5.29p . respectively. The 
latter are in remarkably good agreement with the coefficients obtained from the fit. In addition, 
we give the value -F(z/, 1) which, for all v considered, is compatible with the expected value 
/(I) = 1/2. 



V 


/(oo) 


CO (fit) 


do (exact) 


do (fit) 


F(v,\) 


1.10 


0.548110 


0.548100 


-0.089218 


-0.078176 


0.496148 


1.50 


0.711079 


0.711023 


-0.326980 


-0.319507 


0.488700 


4.20 


1.400983 


1.400000 


-1.049451 


-1.043071 


0.483927 


9.60 


2.418856 


2.413256 


-2.042131 


-2.028276 


0.485305 


20.4 


4.266412 


4.264102 


-3.866293 


-3.858392 


0.489379 
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1/ 










1.10 


-0 564827 


-0 224208 


267706 


066029 


1.50 


-0 745179 


-0 189914 


375405 


231296 


4.20 


-1.613371 


-0.470833 


1.251032 


0.567956 


9.60 


-3.397554 


-3.514617 


3.215461 


3.431349 


20.4 


-5.850448 


-9.119535 


5.615018 


9.110232 



Concerning the values on the last table we see that, for example, the sum cq + ci which we expect 
to be vanishing from the large-n expansion, indeed gives relatively small numbers compared to cq 
and ci, especially for small v. This is consistent with our working precision which is considerably 
reduced for sub-leading terms. A plot of the function -F(n, 1.5) as well as of the numerical values 

of n/^ (0, n) for n = 2,3, ...,70 and v = 1.5 from (j5.7p can be seen in Fig. 2. The difference 

|n/-| (0, n) — F[n, 1.5)| < 10~^ for all n = 2, 3, . . . , 70. In fact, this holds for all finite values of 

considered here. 




Figure 2: The functions F(n, 1.5) and G{n) and the corresponding values of n/-| (0,n) for 

2 < n < 70. 

As we can see from the expansion, the term n^~^/'^ will give a very important contribution 
for v large, so that nf{n) should become linear in n very quickly for v close to its minimum 
value 1 and very slowly for large values of z^. In fact, in the limit ^ cxd, a similar analysis 
as carried out above can be performed (although some complications appear). Both a{j,n) and 
b{j, n) have well-defined v = oo limit which commutes with the 71 = 00 limit, giving 

a(if - ^ hi,? -£ii±Z)! (5 32) 

The leading term at large n, then, can be obtained from 

n/2 n/2 
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where 7 is Euler's number and ^{z) is the logarithmic derivative of Euler's gamma function 
^{z) = T^z)' /T(z). For n large the function ^ (1 + n/2) behaves as 



so that the large-n expansion at = 00 starts with 

lim n/(n) ~ — nlogn 
u^cxD 27r 



(5.34) 



(5.35) 



(where we used b{l/2)'^^^ = S/vr). It is quite interesting to note that this logarithmic behaviour 
can also be obtained from the expansion (j5.30p by simply taking the limit 1^ ^ 00. Both cq and 
do are divergent proportionally to v, and these divergencies cancel out. The expansion in 
of n^~^^'^ then gives the correct logarithmic term. A more precise analysis of cq follows from 



Y.-ij)'bUf ™ ( a(,)^6(j)^ 



+ 



C 1 + - +0(1) 



(5.36) 



where the zeta function contains the linear divergence at large i'. The first term can be evaluated 
exactly: 



E ( «(^')^^(jf 



r(j? 



3 



2 

— ^ lim(2xi^(x) + log(l — x)) 



(5.37) 



(where K{x) is the complete elliptic integral of the first kind). From this, the = 00 limit of 
()5.30p gives the conjecture 

lim 2iTnf{n) = fn + — + + . . .Vog n + (7 + log 32 - log7r)n + — + ^ + . . . (5.38) 

We will now fit the values of nfj^ (0,?!) obtained numerically to a function of the type: 



G{n) = eo(l + n) log(n + 1) + ei log(n + 1) + + + + 1) + A + 



1 + n 



n+l 



(5.39) 



Similarly as -F(z^, ri), the function G'(n) has a large n expansion of the type (j5.38p . The coefficients 
Co and /q from the fit, as well as their exact values from ()5.38p are given in the following table 



eo (exact) 


eo (fit) 


/o (exact) 


/o (fit) 


ei 


62 


/i 


/2 


0.159155 


0.158971 


0.461266 


0.462474 


-0.183370 


-0.124718 


-0.615249 


0.248293 



so that G(l) = 0.4838988. Fig. 2 shows the function G(n) as well as the numerical values 

of n/-| (0, n) from (j5.7p . The agreement between the two sets of values is very good. More 

precisely (0, n) — G{n)\ < 10~^ for all n = 2, . . . , 70. 
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5.1 Discussion of the large-n results 

We notice from the graphs and from the fitting functions that the value of 1/2 at n = 1 is 
perfectly in agreement with smooth and convex analytic continuations agreeing with the given 
values at n = 2, 3, . . ., in all cases analysed. Also, the large-n expansion makes it clear that no 
oscillatory terms appear. We also remark that the leading large-n behaviour is linear in all super- 
renormalisable cases < oo, and is proportional to nlogn in the marginally renormalisable case 
V = oo (this is still a non-conformal case, where a scale appears by dimensional transmutation). 
This seems to suggest that the large-n behaviour is related to the type of perturbing field, 
and we recall that logarithmic functions (of the cutoff) usually appear with naively marginal 
perturbations. Moreover, the full large-n expansion in the finite-z^ case seems to encode some 
characteristics of the perturbing field, as it also has a part with non-integer powers of n (although 
these powers are not in linear relation with the dimension of the field). Also striking and pointing 
towards a perturbation-theory origin of the large-n behaviour is the fact that the limit u ^ oo 
can be taken directly in this large-n expansion, although it is expected to be only an asymptotic 
expansion, and the limit u ^ oo cannot be taken in the various intermediate steps separately. 
We recall that a good regularisation scheme of the sine-Gordon model at = oo is a "non- 
perturbative" dimensional regularisation, whereby u is taken finite (thus changing the dimension 
of the perturbing field) and sent to infinity at the end of the calculation (see for instance |34j ) . 
This leads to logarithmic functions in a similar way to that by which log n appeared here. 

6 Conclusions and outlook 

This paper provides an extension and generalisation of the work [20j to integrable QFTs with 
backscattering. The main conclusion following from [20] and the present work is that the form 
of the leading correction to the bi-partite entropy at large distances of any integrable QFT is 
"universal" in the sense that it does not depend on the precise form of the scattering matrix (in 
particular, whether or not it is diagonal) but only on the mass spectrum of the model. 

The argument leading to this result, developed in [20], uses the well-known replica trick, 
which involves the partition function on Riemann surfaces with two branch points and n sheets, 
related to the correlation function of twist fields. We obtained the first large-distance (distance 
between the branch points) correction to this partition function in the sine-Gordon model. This 
correction presents an interesting large-n asymptotic, which seems to be in relation with the type 
of perturbing field giving the massive theory. It would be very interesting to find a satisfying 
explanation for this large-n behaviour, in particular in the super-renormalisable case, where the 
leading behaviour is just linear. 

From the replica trick, the main subtlety in obtaining the entanglement entropy is the ana- 
lytic continuation in n, the number of Riemann sheets. In the present paper we attempted to 
provide support for the results using the geometrical picture of conical singularities with angle 
27rn instead of branch points, and appealing to some convexity properties as function of n. How- 
ever, a more precise statement about the uniqueness of the analytic continuation we obtained is 
still missing. 

It would be interesting to investigate higher order corrections to the entropy, that is, higher 
particle form factor contributions to the correlation function ()2.14p and, in particular, determine 
whether or not they are of a similar "universal" nature as the leading correction already obtained. 
This would involve the computation of higher particle form factors of the twist fields. Two- 
particle form factors were obtained here by solving a set of form factor consistency equations, 
but they could also be obtained using Lukyanov's method of angular quantisation [35j, as was 
done in the sinh-Gordon model in [20j . This may be quite helpful for higher particle numbers. 
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Also, investigation in other models may help decipher the properties of form factors of twist 
fields, for low-particle numbers as well. The geometrical meaning of the field may provide 
simplifications. For instance, can there be non-zero one-particle form factors? 

Another outstanding point is the generalisation of the approach employed here to the com- 
putation of the entanglement entropy of quantum systems consisting of multiply disconnected 
regions. Such an extension could provide an alternative way to proving certain natural general 
properties of the entanglement entropy shown to be valid in conformal field theory and 
recently discussed in a more general framework in [35]. 

Acknowledgments: We are grateful to John Cardy for discussions. 
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